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Abstract 

A self-consistent tlieory of both spin and charge fluctuations in the Hubbard 
model is presented. It is in quantitative agreement with Monte Carlo data at 
least up to intermediate coupling {U ~ 8t). It includes both short- wavelength 
quantum renormalization effects, and long-wavelength thermal fluctuations 
which can destroy long-range order in two dimensions. This last effect leads to 
a small energy scale, as often observed in high temperature superconductors. 
The theory is conserving, satisfies the Pauli principle and includes three- 
particle correlations necessary to account for the incipient Mott transition. 
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The model introduced initially by Hubbard |l| for itinerant magnets is now widely used 
for High- Temperature SuperConductors (HTSC) and other materials with strong inter- 
electron interaction. Despite the apparent simplicity of the model, its properties remain 
poorly understood in the strong to intermediate coupling regimes relevant for HTSC. Much 
experimental information on the magnetic fluctuations of these materials is now available 
from neutron scattering and nuclear magnetic resonance. A ubiquitous feature of the data is 
the presence of an unexplained small energy scale. The one-band Hubbard model near half- 
filling should contain this feature if it is the correct model for HTSC. Previous explanations 
of the magnetic correlations taking into account short-range quantum correlations ( T- 
matrix effects) explain most of the Monte Carlo data except in the experimentally relevant 
regime. No explanation of the charge structure factor has appeared. 

In this letter, we present a simple self-consistent approach to the two-dimensional Hub- 
bard model that gives, without adjustable parameter, quantitative agreement with Monte 
Carlo data for spin and charge structure factors and susceptibilities at all fillings up to quite 
strong coupling. The approach takes into account not only the short-range quantum effects, 
but also the long-range thermal fluctuations that destroy antiferromagnetic long-range or- 
der in two-dimensions at any finite temperature (Mermin- Wagner theorem). This is the 
key physical ingredient which leads to a small energy scale, and associated large correlation 
length, in the magnetic fluctuations. Previous approaches which included the long-range 
thermal fluctuations in the 2D Hubbard model were never applied to the incommensurate 
case relevant for HTSC. Furthermore, they are based on mode-coupling theory which 
neglects charge fluctuations and does not include the effect of short-range quantum correla- 
tions, which are important not only for a quantitative description of the model but also for 
determining the nature of the ground state. 

We flrst present our approach, discuss physical consequences, and flnally compare with 
Monte Carlo data. We consider the one-band Hubbard model with on-site repulsive U . Our 
approach is motivated by the Local Field Approximation (LFA), which was successful in 
the electron gas P|. We start from the equation of motion for the particle-hole operator 
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Po- (j., l'^ = ct^cp in a weak external field (j)^^ which is coupled to the partial density 
operator nj^^ = p„ (^)^ — 0- After simple transformations, the term that contains 

the interaction in the equation of motion for p„ U,^') is of the form, 



U [{p.{lJ')nr^_J-{p^{ir')nr,^_J]. (1) 

All operators contain the same time label, which is not explicitly written. The usual Ran- 
dom Phase Approximation (RPA) corresponds to the neglect of two-particle correlations, 
namely one approximates {pa (l, /') P-o- ^ ) by (po- (j,, l'^ ) (p_o- (^7 ^ )• This is clearly a poor 
approximation for on-site interactions because, as can be seen from ([^), three of the four 
creation or annihilation operators in the correlator {p„ (^'^') P-^ (j^^) same lat- 

tice site /. There is thus a strong correlation between two particle-hole pairs even when 
| / — /'I >> 1 (the lattice constant is taken to be unity). We make use of this specific feature 
of the on-site interaction and neglect the dependence of the correlation coefficient on the 
lattice index which appears only once in the two-particle correlator. Mathematically, our 
ansatz is, 

(p. (r, I') nr^_J = gn {l, 1} (p. {i, f) ) {nr^_J (2) 

where 

— (Mwj) — 

is the pair correlation function between electrons of spin a and a' on the respective lattice 
sites / and 

It is important to realize that the pair correlation function in (|l|) and cannot yet be 
taken equal to its equilibrium value (l,l;tj 7^ g^i (0) because of the weak external field 
(j)^^. In linear response the most general form for g-^i is: 



gn {l, I; t) = gn (0) + E / dt' sll^^\t>^) ^^^r' (t')) 
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where nr{t) = rif^ (t) + nr. (t) and we use that in the paramagnetic state ''^ti 



5^{n'-}'^ie\ = 5'n {} ~ ^' ~ Note that the terms describing the response of g^i (l,l;t 
to the external field enter the equation of motion only in a form that is symmetric in spin 
indices. Since the 2;— component of spin = — ni is antisymmetric in rig- and all equa- 
tions are linear, it immediately follows that (7'^ (j, — l',t — t'^ enters the equation for charge 
but not for spin. This important simplification occurs because the Pauli principle precludes 
the appearance of terms like 7^ in the case of on-site interaction. After standard 
transformations (see for example 0) the spin and charge susceptibilities have the RPA form 
but with a renormalized effective interaction, which is different for spin (Usp) and charge 

Usp = 9n (0) U- t/e, = {gu (0) + 6g^^ {u, q)) U (4) 

where 5g^i (u;, q) = g'^^ (cj, g) |, g'^^ {uj, q) is Fourier transform of g'^^ (j. — l',t — t'^ and n is the 
band filling (half-filled case corresponds to n = 1). 6g^i {'^i<i) is a three-particle correlation 
function, so that further simplification will be needed to calculate the charge susceptibility 
Xch iyJ, q)- However, no further approximation is needed for the spin susceptibihty Xsp iyJi 
Indeed, due to the Pauli principle, g^^ (0) = so the spin part of the problem may be closed 
by using g^iiS^) = —2gsp{0) with gsp{l,l') = ((7||(Z,Z') — g^i{l,l'))/2 and the Fluctuation- 
Dissipation Theorem (FDT) for spin: 

^?«p(M') = ^/(i?.[^sp(g-)-i]e-M 

^sp „ ^iujm l-{Usp/2)xo 

The first equation is the definition of the spin structure factor; the second is a convenient 
form of the FDT with temperature T and bosonic Matsubara frequencies ium- The integral 
is over the first Brillouin zone. The definition of xo {'i^m, q) is the same as in Ref. [0. 

For 5g^i {uj,q), we use the simplest ansatz, namely, that it is a constant 6g^i, which we 
determine self-consistently using the Pauli principle g^^ (0) = 0, the definition for the static 
charge structure factor 



9cH (/,/') = 1 + - / 7^ [ScH {q) - 1] e^^^('-'') (6) 

and the FDT for charge. We thus have a simple theory with only two parameters g^y (0) 
and 5g^i that are found self-consistently. It can be explicitly checked that charge, spin, and 
energy are conserved. As with all self-consistent theories, the usefulness of the approach can 
be judged only a posteriori by comparison with numerical or exact results. We provide such 
comparisons later in the paper. 

The absence of a magnetic phase transition at any finite temperature in 2D follows 
immediately from the above approach. Define the mean-field critical value of U by f/m/,c = 
2/xo (0, ^max) where xo (i^m, ^nax) is the susceptibility of non-interacting electrons and gmax 
is the value of q at which the static susceptibility has its maximum. If Ugp = g-ii{0)U in 
(|) was large enough for the transition to occur, namely 6U = \Umf,c — fi'ti (0) ^] — 0; 
then the g— integral for the static susceptibility {um = 0) in the expression for gspiS^) (0) 
would diverge logarithmically so that g^i (0) = —2gsp{^) would become negative, in obvious 
contradiction with 5U = 0. Hence, in our approach, magnetic fluctuations always push the 
value of g-^i (0) away from its critical value g^^^ (0) = Umf,c/U at any finite temperature. 
Furthermore, for a wide range of values of U > Umf,c and T < Tm/,c , the system will be 
quite close to magnetic instability {6U ~ 0), providing the basis for a generic explanation of 
the small energy scale observed in HTSC. In the regime in which the temperature is larger 
than this small energy scale, the correlation length grows exponentially ^ oc oc e'^""^*/-^ 
(^^^ = SU/Umf,c), reflecting the logarithmic divergence mentioned above. This is typical 
behavior for systems at their lower critical dimension. When a real quasi-two-dimensional 
system enters this regime, small three-dimensional effects can easily stabilize a long-range 
order. 

We digress briefly to speculate on how the phase diagram of a system with weak three- 
dimensional effects would then look. We neglect effects, such as disorder, which may be- 
come important when the small energy scale SU appears. We define a quasi-critical tem- 
perature Tqc as the temperature at which the enhancement of the magnetic susceptibility 

5 



= Xsp (0, ^nax) /Xo (0, i^nax) IS of Order 500. Calculations for the nearest-neighbor Hubbard 
model show that increases by an order of magnitude when the temperature is reduced 
below Tgc by as little as Tg^ — T ~ 0.01 (all energies are in units of the hopping integral t). 
In our theory, the emerging long-range order is determined by the position of the maximum 
of Xo (0, ^. At T = 0, as soon as n 7^ 1, two-dimensional Fermi-surface effects ^ lead to 
a maximum with a cusp-like singularity in the (tt, tt) — (tt, 0) direction. The situation is 
different at finite temperature where the maximum can be at g^ax = (tt, tt) even for n < 1. 
Fig. 1 shows a rough magnetic phase diagram obtained by approximating Tc by Tgc (with 
U = 2.5). The inserts show the dependence of gj,,max (^nax = (ti"? %,ma.x) ) and of the enhance- 
ment factor on temperature for three different fillings. The filling n = 0.93 is marginal in 
the sense that the shift of g^ax from {it, tt) occurs when the enhancement is already quite 
large. We would expect then that long-range order would be antiferromagnetic (AFM) for 
n > 0.93 and incommensurate spin-density wave for 0.89 < n < 0.93. 

Let us now discuss the limit T — > 0. In Fig. 2, we show the renormalized spin interaction 
Usp as function of the filling n for U = 2.5 and U — * 00, together with Umf,c (n). Remark- 
ably, at low filling the value of Ugp saturates to Ugp ~ 3.2 as the bare interaction increases 
U —>■ 00. This quantum effect was anticipated by Kanamori [0], ^ who argued that the 
largest value of Usp is proportional to the kinetic energy cost to put a node in the two-body 
wave-function where two electrons overlap. At a sufficiently large filling, when Usp{n) starts 
to follow Umf,c (n) for T — 0, the paramagnetic state has an instability at exactly T = 0. 
The existence of an upper limit for Usp leads to the existence of a lower limit for the filling 
nc,min ~ 0.685 below which there is no magnetic phase transition at any U. This, in turn, 
means that only spin-density waves with q = (vr,gy), qy/n C [0.74,1] are possible. In 
particular, the ferromagnetic state does not exist in the Hubbard model on a square lattice 
and there is also a temperature T 0.5 above which there is no exponential regime for 
^{T) at any U where our theory applies. 

In the rest of this letter, we compare our theoretical results for infinite lattice with Monte 
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Carlo simulations. Finite-size effects are expected to become important in simulations when 
the correlation length becomes comparable with system size . The largest size effects are 
thus expected in Sgp (tt, tt) at half-filling (n = 1). This is shown in Fig. 3. The Monte Carlo 
data follow our theoretical curve (solid line) until they saturate to a size-dependent value. 
We checked that finite-size effects for Sgp (q) away from the antiferromagnetic wave vector 
q = {it, tt) are much smaller, so that, even for 8x8 systems, theory and simulations agree 
very well for all other values of q (not shown). Finite-size effects in the spin structure factor 
Ssp (q) are not too important away from half-filling for the parameters shown on Fig. 4a-d. 
Obviously though, finite-size simulations cannot reproduce the small incommensurability 
captured by our theory at n = 0.8 on Fig. 4d. 

Figures. 4b and 4c show that, even for relatively strong coupling {U ~ 8), the theory 
agrees very well with both spin Sgp (g) and charge Sch (q) structure factors. However the 
theory should eventually break down for [/ — oo. This can be seen in the half- filled case 
from the fact that, for f/ ^ oo, the antiferromagnetic susceptibility remains constant in 
our theory while mapping to the Heisenberg model with J = At'^/U shows that it should 
decrease with U. It seems, however, that large-f/ asymptotic is reached for values of U much 
larger than the bandwidth. 

Fig. 4a shows that our theory reproduces the important qualitative fact that the charge 
structure factor Sch [q] depends on filling in a non-monotonous manner. The decrease 
of Sch (q) towards half-filling is a signature of the incipient Mott transition. The effect can 
be seen because our approach takes into account both three-particle correlations and the 
Pauli principle. Writing the Pauh principle as a sum-rule Sg-[S'c/j (g) + Sgp (q)] = 2T,^So (q), 
the parameter Squ, which partially takes into account three-particle correlations, must in- 
crease close to half-filling in order to reduce Sch (q) and compensate for the increase in the 
contribution of the spin structure factor. 

Our theory also explains the good fit of the dynamical spin susceptibility Xsp {i^^m, q) 
obtained by Bulut et al. using RPA with Uren = 2. Indeed, for f/ = 4, n ~ 0.87 on 
8x8 clusters, our calculations give Usp = 2.05 with very little dependence on temperature. 
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Bulut et al. have also shown that their Monte Carlo data for the self-energy S {iujm, q) can 
be reasonably well fitted by the Berk-Schrieffer formula |Ty] with the same Uren = 2. In our 
approach, the expression for S {ium-, <l) in terms of susceptibilities should come at the next 
level of approximation. Bulut et al. have also fitted a number of experiments in HTSC 
by fine-tuning the value of f/^en close to a magnetic instability {5U ~ 0). In our approach, 
a wide range of bare values of U naturally renormalizes to such a situation. 

In conclusion, imposing the Pauli principle as well as self-consistency through the 
fluctuation-dissipation theorem, we have formulated a simple theory that also satisfies con- 
servation laws and gives, without adjustable parameter, a quantitative explanation of Monte 
Carlo data for both spin and charge structure factors as well as susceptibilities up to interme- 
diate coupling. Both short-wavelength quantum renormalization effects and long-wavelength 
thermal fluctuation effects, which destroy long-range order in two- dimensions, are accounted 
for. The latter effect naturally leads to a small energy scale for a wide range of parameters, 
possibly giving a microscopic origin for the small energy scale observed in experiments on 
high-temperature superconductors. 
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FIGURES 



Fig. 1. Approximate phase diagram of the quasi-2D Hubbard model with Tc approxi- 
mated by Tqc. The insert shows the temperature dependence of qmax {T) and the enhancement 
factor f2 = Xsp{0, qmax)/XoiO, qmax) for three different filhngs n. 

Fig. 2. Filhng dependence of Usp and Umf,c as T ^ 0. 

Fig. 3. Temperature dependence of Ssp^ir, vr) at half-filhng n = 1. The sohd hne is our 
theory; symbols are Monte Carlo data from Ref. 

Fig. 4. Wave vector (q) dependence of the spin and charge structure factors for different 
sets of parameters. Solid lines are our theory; symbols are our Monte Carlo data. Monte 
Carlo data for n = 1 and U = 8 are for 6x6 clusters and T = 0.5; all other data are for 
8x8 clusters and T = 0.2. Error bars are shown only when significant. 
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